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Abstract 

In our recent publication pp we presented an exponential series ap- 
proximation suitable for highly accurate computation of the complex 
error function in a rapid algorithm. In this Short Communication we 
describe how a simplified representation of the proposed complex error 
function approximation makes possible further algorithmic optimiza- 
tion resulting in a considerable computational acceleration without 
compromise on accuracy. 
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Consider the exponential approximation of the complex error function, 
obtained by Fourier expansion method [1] 
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where z = x + iy is a complex argument and 
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is a set of the Fourier expansion coefficients. 

Let us take a term corresponding to index n = apart from the summa- 
tion sign in approximation series ([1]). This leads to the following representa- 
tion 
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The expression inside parentheses of the approximation ([2]) can be rearranged 
as 
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Since (e in * - e~ m7r ) /2i = sin(mr) = and (e m7r + e~ m7r ) /2 = cos (rwr) ee (-1)", we 
can rewrite it in form 
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Ultimately, applying this relation in approximation (T5]) yields 
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Comparing the approximation representations (PQ) and ([3]) with each other, 
we can see that the second one is substantially simplified and we can take 
advantage of it for computational optimization. 

In order to describe shortly and conveniently how the simplified represen- 
tation of the complex error function ([3]) can be appropriately implemented 
in a rapid algorithm, we will use the standard Matlab notations for element- 
by-element array multiplication " . *" and division " ./" with regular and bold 
letters for scalars and arrays, respectively. Thus the algorithmic implemen- 
tation based on approximation fl3]) can be symbolically conformed to 
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where A := r m z, B := e lA , and C := A. * A are the only three arrays 
involved in computation (after initialization of A, the array z is never used 
anymore and can be immediately cleaned from the computer memory). Since 
none of these arrays are dependent on index n, each of them are computed 
just a single time. As follows from this scheme, only the array B requires 
the exponentiation. All other operations are trivial such as summation, mul- 
tiplication and division. Since the calculation of the array B is made only 
once, it has a very minor contribution (at r m = 12 and N = 23 it is less than 
6%) for the total computation time. This explains why the calculation of the 
Voigt /complex error function is de facto as rapid as a rational approximation. 

The notation e lA implies element-wise exponential operation over the ar- 
ray A that returns the array B with same dimension; such element-wise 
functional operations are typical in array programming languages like Mat- 
lab/GNU Octave or modern FORTRAN. 

It should be noted that if the high-accuracy is not needed, the acceleration 
can be additionally gained by reducing the values of r m and N, say to 9 and 
12, respectively. The approximation of the exponential function e lTmZ by a 
rational function is not required as its contribution to the total computation 
time remains insignificant even in this case. 

The advantage of the approximation ([3]) implementation becomes partic- 
ularly evident at an extended size of the input array. Specifically, the com- 
putational test reveals that with size of input array z exceeding 30 million 
elements, the algorithm based on approximation is more than twice faster 
than the one based on approximation (JT|) and, compared to the Weideman's 
algorithm (we used Matlab code in page 1511, Ref. [2]), it becomes faster by 
a factor greater than 5 (see also Ref. p] for comparison). 

Similar to approximation ([T]) the numbers generated by approximation 
02]) match up to the last decimal digits with those generated by well-known 
Algorithm 680 [HE]. This signifies that due to simplified form of the approx- 
imation ([3]) the considerable acceleration in computation is achieved without 
compromise on accuracy. 

All calculations were performed with help of Matlab 7.9.0 on a typical 
desktop computer Intel(R) Quad CPU with RAM 4.00 GB. 
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